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Abstract 

The ionized disc reflection model of Ross & Fabian is 
applied to a sample of five Seyfert 1 galaxies observed by 
BcppoSAX. Good fits were obtained for all the objects in 
the sample. NGC 4051 and Mkn 335 (usually classified 
as NLSy Is) show evidence of a highly ionized accretion 
disc (£ ~ 350 erg cm s _1 and £ ~ 800 erg cm s _1 ). While 
NGC 3783 seems to be a "pure Sy 1" characterized by a 
highly ionized accretion disc (£ ~ 500 erg cm s" 1 ). All 
the sources are consistent with having a face-on disc. The 
accurate determination of the continuum emission allows 
us to investigate the spectral variability in the case of Mkn 
509. 
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1. Introduction 

Evidence for ionized accretion discs in Sy Is was found in 
rece nt observations by X MM and XMM-BeppoSAX (Mkn 

NCG 



509, Pounds et al. 2001 



Mkn 205, Reeves et al. 2001 



5506, Matt et al. 2001 ). This conclusion was reached by 
considering the energy of the Fe Ka line. The ionized disc 
also affects t he whole continuum Ross fc Fabian 1993 , 



hereafter RF, Ballantyne et al. 2001b ). The wide energy 
range of BcppoSAX (0.1-200 keV) and its moderate spec- 
tral resolution (AE/E 8% at 6 keV) allow us to mea- 
sure the intrinsic continuum with good accuracy and to 
simultaneously observe both the Fe emission line and the 
Compton reflection hump. We present preliminary results 
of fitting BeppoS AX spectra of a sample of Seyfert 1 galax- 
ies with the ionized disc model of RF. 

2. The selected sources 

Details of of BeppoSAX observations and the sources are 
shown in Table [l]. The sample includes observations of 
sources with F(2-10 keV) greater than 10~ n erg cm~ 2 
s _1 and exposure time greater than 50 ks (with the excep- 
tion of NGC 4051). This allowed us to obtain a good signal 
in the PDS band. We search all the sources for spectral 
variability. We accepted a source as variable if the proba- 
bility P(x 2 ) that either the MECS (5-10 keV)/MECS (3-5 



keV) or PDS (13-200 keV)/MECS(5-10 keV) hardness ra- 
tio was constant was < 0.01. This allowed us to merge 
together the four short observations of Mkn 509 in 2000. 
The observation of this source in 1998 had a flux higher 
by a factor of 2 from the observations performed in 2000, 
so it is analysed separately. NGC 3783 also showed mod- 
est spectral variability, but the strength of the variations 
in that source is n ot relevant to the present work (see 
Do Rosa et al. 2002] for details about spectral variability). 
Therefore, we analysed the averaged spectrum extracted 
by the total observation. 

3. The Ionized disc model 

The ionized disc model we used to fit the spectra is de- 
scribed in RF and Ross et al. (1999)| . The most important 
quantity in determining the shape of the reflected contin- 
uum is the ionization parameter £ = 4irF x /riH, where F^ 
is the X-ray flux (between 0.01-100 keV) illuminating a 
slab of gas with solar abundances and constant hydrogen 
number density n#=10 15 cm" 3 . The incident flux is as- 
sumed be a power law with spectral photon index T and 
a high energy cut-off at 100 keV. The computed reflected 
spectrum is multiplied for a factor R (reflected fraction) 
and added to the primary continuum. The Fe Ka emission 
line is also included in the model. Larger values of £ indi- 
cate a more ionized disc a nd this will affect t he strength 



and width of the Fe line flMatt et al. 1993| , Matt et al 



1996 ), and the absorption edges. The soft emission from 
the accretion disc is not included in the models. 



4. Spectral fitting 

Each source spectrum was fitted with the RF model de- 
scribed in the preview section (the best fit parameters are 
listed in Table |). We fitted LECS (0.15-3 keV), MECS 
(1.5-10 keV) and PDS (13-200 keV) data simultaneously. 
The cold absorption is kept fixed to the galactic value. 
Rclativistic smearing was also taken into account by blur- 
ring the model with a kernel of Schwarzschild metric. This 
introduces the following fitting parameters: Ti n and r ou t- 
inner/outer emitting radius of the annulus, (3: radial emis- 
sivity parameter of the disc and i: inclination angle. We fit- 
ted all the spectra also with a standard cold disc reflection 



model (PEXRAV in XSPEC, |Magdziarz fc Zdziarski 1995 . 
The best fit parameters listed in Table |3|) in order to: (a) 



Proc. Symposium 'New Visions of the X-ray Universe in the XMM-Newton and Chandra Era' 26-30 November 2001, 



2 



Table 1. The Sample 



Source 


Date 


Redshift (z) 


N Gal 

(in 10 20 cm" 2 ) 


cts (2-10 keV) 
(s- 1 ) 


(MECS) 


T^obs 

r (2-10fceV) 

(erg cm 2 s" 1 ) 


Ref 


Mkn 335 


1998-Dec-10 


0.025 


3.7 


0.090 ± 0.001 


86637 


0.8X10" 11 


3ianchi et al. 2001 


Mkn 841 


1999-Jul-30 


0.036 


2.2 


0.130 ± 0.001 


87793 


1.3xl0 _li 


3ianchi et al. 2001 


NGC 3783 


1998-Jun-6 


0.009 


9.6 


0.663 ± 0.001 


153870 


6.0xl0~ n 


3e Rosa ct al. 2002 


NGC 4051 


1998-Jun-28 


0.002 


1.3 


0.203 ± 0.004 


11616 


1.9xl0~ n 


This proceeding 


Mkn 509 


1998-May/Oct 


0.034 


4.4 


0.603 ± 0.004 


87762 


5.6xl0" n 


Perola et al. 200C| 


Mkn 509 


2000-Nov-3/24 


0.034 


4.4 


0.278 ± 0.002 


152380 


2.6xl0~ ri 


This proceeding 



Table 2. Results I: The ionized disc reflection model. The disc emissivity law is e(r) oc r' 3 , with (3 =-3. The inner radius 
is kept fixed: ri n =6 r g . with r g — Gm/c 2 the gravitational radius. 



Source T log (£) R r otlt (r 9 ) i(degrees) X 2 /dof 



*Mkn 335 


2.08±££ 


2 ' 9± ii : L 


9 i , U.V 
^•■•■±0.5 


*6.9 


29.5±£i 


97.8/94 


Mkn 841 


i.97±8:8J 


1.38± ; 28 




>40 


<20 


115/109 


*NGC 3783 


1.70±0.01 


9 7 Q 1 0.04 


0.49±H! 


*6.9 


35.6±?; 2 


104/101 


*NGC 4051 


i.69±8:8S 


2.65± ; g 


1.18±g;| 


7±i 


23.5 ii7 5 2 


99.5/106 


Mkn 509 (1998) 


1.90±0.01 


1.36±8;it 


0.83± ;22 


34±26 


32±10 


120/143 


Mkn 509 (2000) 


1.72±0.02 


1.43±0.19 


1.44±°;i 


*10 


ll±f 


128/143 



• For these sources the parameters are the best fit model has to take into account the complex soft X-rays spectra. A warm 
absorber component is therefore added to the RF model (see discussion in Section . 

* Parameter fixed at that value in the fit. 



separately determine the Fe Ka parameters and (b) check 
the relative importance to employ the ionized disc model. 
All the SAX data were fitted using the XSPEC package 
11.1. All the quoted uncertainties correspond to 90% con- 
fidence interval for one interesting parameter (Ax 2 =2.71). 
Each model tested was multiplied by a normalization con- 
stant in order to take into account possible miscalibra- 
tions between the different instruments. We allowed the 
PDS/MECS normalization to vary between 0.77 and 0.95, 
while the LECS/MECS normalization ratio was running 
between 0.7 and 1 ( Fiore, Guainazzi fe Grandi 199S ). 



5. Results: summary and conclusions 

The RF ionized disc reflection model provides an ad- 
equate fit for all observations examined. When a cold 
disc reflection model in employed to fit the spectra 
the general trend is a higher value of the spectral in- 
dex and, in addition, a soft X-ray emission component 
(due to the disc emissivity) in required. 
In the case of NGC 3783, Mkn 335 and NGC 4051 the 
complex low energy spectra (characterized by warm 
absorber(s) and strong soft excess) require a more de- 
tailed study. 

The complex structure of the warm gas in NGC 3783 



was investi gated in the Chandra observation (Kaspi 



lg(iVff) = 22.1, and differing each other an order of 
magnitude: lg(Z7) = -0.745, and lg({7) = -1.745. 
With the BeppoSAX-LECS energy resolution the warm 
absorber is fitted through the detection of the deep 
OVII and OVIII K-shell absorption edges produced in 
the low ionization gas. We added to the ionized disc 
reflection model in NGC 3783 two absorption edges 
at energies 0.74 keV (r=0.71) and 0.87 keV (t=0.53) 
respectively, coincident with the OVII, OVIII K-shell 
edges energy, and a third absorption feature at 1 keV 
probably due to a blend of many iron L-shell absorp- 
tion lines originated in the second more ionized gas 
(Kaspi et al. 2001). In the best fit model we included 
also a narrow Fe line with the intensity fixed to the 
value observed by Chandra (Ijr e _ narro w =6.6 x 10 -5 
photons cm" 2 s" 1 , Kaspi et al. 2001). As shown in 



Table |2| the fit for the ionized disc reflection model is 
good (x 2 /dof=104.2/101). 

In the case of Mkn 335 the best fit model in the whole 
BeppoSAX energy range includes a warm absorber 
component characterized by OVII and OVIII absorp- 
tion edges (r=0.67 and t=0.53 respectively). The strong 
soft excess observed in the low energy spectrum can be 
accounted by the emission from a highly ionized disc 
(lg(£)= 2.9, see Table §). 

NGC 4051 has a warm absorber well known for its 



et al. 2001 ). The model consisted of two absorption- 
emission components with the same column density 



spectral complexity ( Guainazzi ct al. 1996 ). Neverthe- 
less the BeppoSAX low energy spectrum was well fitted 
adding to the RF model an absorption edge at energy 
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Table 3. Results II: The iron line. 



Source 


r 


(KeV) 


a or T out 
(keV) or (r 9 ) 


EW 
(eV) 


X 2 /dof 


OMkn 335 


2.04±8:& 


*6.4 


6.1±g 


660±150 


126/123 


Mkn 841 


9 19 , 0.10 


*6.4 


g.o±i:§ 


260±|g° 


120.8/119 


NGC 3783 


i.87±8:8g 


6.39±0.09 


o.36±8:81 


190±47 


100.7/102 


NGC 4051 


1.65±0.10 


6.2±?;I 


<i 


135±200 


102.3/106 


Mkn 509 (1998) 


1.97±0.03 


6.88±0.01 


*0.4 


115±65 


101.3/141 


Mkn 509 (2000) 


1.75±0.06 


6.57±8:l! 


0.55± f 2 ,;? 8 


149±100 


106/126 



* Parameter fixed at that value in the fit. 

^ The Fe line profile is modelled adding to the ralativistic disc line a second narrow gaussian component at E=6.4 keV. 



0.88 keV (t=0.74). Also in this case the strong soft 
excess can be accounted by the emission from a highly 
ionized disc (lg(£)= 2.65, see Table |f). 
We stress here that Mkn 335 and NGC 4051 are clas- 
sified as a Narrow Line Seyfert 1 galaxies, and in this 
class of sources high ion ized disc are often observed 
flBallantync ct al. 2001a] ) . 

The RF model a pplied on a sample of NLSy Is ob- 
served by ASCA (Ballantyne et al. 2001a), shows that 
these sources might posses an ionized disc. Most of 
the Sy Is we analysed observed by BeppoSAX show 
evidence of an ionized disc. The ionization parameters 
are anyway much less that that observed in NLSy Is, 
not surprisingly as in latter class of objects are com- 
monly believed to accrete at high rates. NGC 3783 is 
the "pure Sy 1" in our sample that seems to be char- 
acterize from a highly ionized accretion disc (lg(£) = 
2.73). 

All the sources well reproduced by a RF model are 
characterized by a face-on disc with i < 40° (consistent 
with the ASCA observations of Sy Is, Nandra et al. 
1996) and with radial emissivity parameter /3=-3. 
The flux variation of a factor of 2 detected between the 
two observations of Mkn 509 (taken 2 years apart), can 
be completely explained with a change of the intrinsic 
slope Ar ~ 0.2. This behaviour is fully consistent with 



that observed in other Sy Is (NG C 5548, Nicastro ct al. 
2000; NGC 4151, |Piro ct al. 2002|; IC 43 29a, [Done ct al. 
200 1 NGC 7469 iNandra et al. 20001 MCG-6-30-15 



Edelson 2001 



Vaughan 

200g), and with the F , 



Dc Rosa et al 



NGC 3783 
vs a relationship expected in 
a Comptonization model for the intrinsic continuum 



emission (Petrucci et al. 2000 and references therein) 



We are extending this analysis to a wider sample (20 
objects) of sources observed by BeppoSAX. 
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Figure 1. In each plot we show LECS (black), MECS (red) and PDS (green) data (upper panels) and data/model 
ratio (lower panels) for the ionized disc reflection model (the best fit parameters are listed in Table ^). In the case 
of NGC 3783, Mkn 335 and NGC 4051 the spectral analysis must take into account the complex soft X-rays spectra 
characterized by warm absorber(s). These features are included in the best fit model (see discussion in Section^). The 
soft excess observed at E < 2 keV (particularly strong in the two NLSy Is NGC 4051 and Mkn 335), can be accounted 
by the disc emissivity. 



